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Abstract 

A numerical model of spontaneous decay continuously monitored by a distant de- 
tector of emitted particles is constructed. It is shown that there is no quantum Zeno 
effect in such quantum measurement if the interaction between emitted particle and 
detector is short-range and the mass of emitted particle is not zero. 
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1 Introduction 

2 ■ 

The Quantum Zeno Paradox (QZP) is a proposition that evolution of a quantum system is 
stopped if the state of system is continuously measured by a macroscopic device to check 
whether the system is still in its initial state HJ. QZP is a consequence of formal ap- 
plication of von Neumann's projection postulate to represent a continuous measurement 
as a sequence of infinitely frequent instantaneous collapses of system's wave function. It 
was shown theoretically || and experimentally || that sufficiently frequent discrete active 
measurements of system's state really inhibit quantum evolution. This phenomenon was 
named 'Quantum Zeno Effect' (QZE). But the question about possibility of QZE during 
true continuous observations is not quite clear up to now. 

A true continuous measurement of quantum system's state takes place during observation 
of spontaneous decay by distant detector of emitted particles (another examples of continuous 
measurement of decay are presented in papers || ||, [7], §, [|). Let us consider a metastable 
exited atom surrounded for detectors to register an emitted photon (or electron) when the 
exited state of atom decays to the ground state. While the detectors are not discharged, the 
information that the atom is in its exited state is being obtained permanently, therefore the 
system's state is being measured continuously. Could the presence of detectors influence on 
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the decay constant of exited atom? If so, this would be the QZE in true continuous passive 
measurement. 

It is impossible to describe this kind of continuous measurement by a sequence of discrete 
wave function collapses as was proposed in seminal works [[3], §]. Such approach leads to 
the explicit quantum Zeno paradox, not effect. Instead, a dynamical description of such 
measurements was elaborated in the number of works [|l^, 11, 12, [l3|. In this approach object 



system (atom), radiation field (or emitted particle), and device (detector of particle) are 
considered as subsystems of one compound quantum system. The results of papers |10|, pTT| 
were mainly qualitative. The explicit expression for decay constant perturbed by given 



interaction W of emitted particles with detector was obtained in [|T|, [13| . This expression is 

r = 2tt J duM(uo)A(uj - S ). 



1) 



In Eq. ([!]) M(oj) is the sum of all transition matrix element squares related to the same energy 
of emitted particle u; So is the expectation value of final energy of emitted particle. The 
function A(lj — So) describes the influence of observation on the decay constant. Without 
detector, i. e. W = 0, the function A(u> — So) transforms to Dirac's delta-function 5(uj — So) 
and Eq. ([[J) transforms to the Golden Rule |L2], |13| . It was supposed that the meaning of 



A(u — So) is an energy spreading of the final states of decay due to time-energy uncertainty 
relation and finite time-life of emitted particle until scattering on the detector^]. Then one 
can suggest that observation influences on decay in accordance with the following sequence: 
The faster detector, the shorter emitted particle time-life, the wider A(uj — So), the stronger 
perturbation of decay constant. 

It is clear that strong interactions W is needed to obtain QZE. Hence, W is essentially 
nonperturbative in this problem. This feature determines the main difficulty of calculations 
of function A(u> — S ) in Eq. (P and, consequently, the perturbed value of decay constant. 
Particularly, in paper [T2|] we supposed that QZE explains strong inhibition of 76eV-nuclear 
uranium-235 isomer decay in matrix of silver ||14|| . However, we had to restrict the consider- 
ation only by a qualitative analysis of Eq. (0) for this case because of difficulties of function 
A (a; — S ) calculations. 

Since it is difficult to study realistic physical systems, it is reasonable to start with 
some simplified models to calculate the function A. The aim of the present paper is strict 
and complete numerical investigation of Eq. (TT]) for a simple but not oversimplified model 
system. We derive Eq. ([]]), then introduce one-dimensional three-particle model of continuous 
observation of decay, then describe the numerical computation scheme for this model, and 
finally discuss results of calculations. 



2 General considerations 

In this section a derivation of Eq. (H) and other formulae to construct our numerical model 
are presented. Derivation of Eq. ([!]) is simplified in comparison with our previous papers 

E3E3- 

1 This supposition was confirmed in a case of problem of decay onto an unstable atomic state This 
problem is close to problem of observation of decay by distant detector. 
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Figure 1: The compound system S = X ®Y <g> Z is the model of continuous observation of 
exited state \X e ) decay. 

Let a compound system S = X ® Y <g> Z consists of three subsystems X, Y, and Z 
(Fig. 0)0. The system X ("atom") decays spontaneously from the initial exited state \X e ) 
to the ground state \X g ) emitting a particle Y ("electron") due to interaction V between 
systems X and Y . This process is similar to autoionization decay of excited atomic state, 
but it is possible to suppose another nature of systems and interactions. The particle Y 
is initially at the ground state \Y ) (electron is on the bounded state in atom) and then 
transits to continuum \Y(r], Ey))- Here Ey is the energy of the state in the continuum and 
r) represents all other quantum numbers. Particle Y inelastically scatters on third system 
Z ("distant detector") due to interaction W between Y and Z. As a result, the system Z 
transits from the initial ground state \Z ) to the continuum \Z((, Ez))- This transition is 
considered to be a registration of decay. We consider that the interaction V does not effect 
on system Z, the interaction W does not effect on system X and the systems Y and Z don't 
interact in their ground states. Therefore, we have 

V = V X Y®Iz] W = Ix®W YZ ] W yz \Y Z ) = 0, (2) 

where Ix and Iz are unit operators in the Hilbert spaces of corresponding systems. The 
Hamiltonian of whole system is 

H = H + V + W, (3) 

where 

H = H x ® I Y z + H Y ® I XZ + H Z ® Iyx 

with obvious notations. 

2 S = X®Y ®Z means that the Hilbert space of system S is a direct product of spaces of systems X,Y, Z. 
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The initial state of system S at the initial moment of time T = is 

|#0> = \X e ) ® \Y ) ® |^o) = |X e F ^0). 

Let us introduce the first order correction to the eigenenergy of state |^q) due to interaction 
V: 

5V = (V \V\V ) 

and renormalized unperturbed Hamiltonian Hq and renormalized interaction V: 

H' = H + SV \* )(* \; V' = V- <^ |*o)(*o|- 
Then the Hamiltonian Eq. @ may be rewritten as 

H = H'q + V + W. 

The initial state |^o) is an eigenstate of the Hamiltonian H' with the eigenenergy 

£q = E e x + Ey + E% + 6V . 

The interaction V is considered to be a small perturbation, but interaction W is not 
small. To obtain the decay constant of the exited state \X e ) it is necessary to solve the 
Shrodinger equation for the whole system S. It is impossible to construct the perturbation 
theory for W, but it is possible for V. Therefore, let us introduce the interaction picture as 
(ft = l): 

|tf 7 (T)) = e^ H 'o +w ^ T \^(T)}, |^(0)) = |tf ) (4) 
VUT) = e i{H o +w)T V'e~ iiH o +w)T . 

Then the Shrodinger equation reads as 

|*/CT) = |* > ~i [ T V{(t)\Vi(-t))dt (5) 
Jo 

The solution of Eq. (|5|) in the second order of perturbation theory with respect to V is 

|tf/(T)> = |* )-i [ T V;(t)\V }dt- C dh r dt 2 Vi(ti)Vi(t 2 )\^Q). (6) 
Jo Jo Jo 

Let F(T) be no-decay amplitude 

F(T) =e< T (m Q \m(T)). 
It follows from Eq. (g) and Eq. (§) that 

F(T) = 1 - C dt, [ h rft 2 <*o I Vr(ti) V?(t 2 ) | * > • (7) 
Jo Jo 

For the initial region of exponential decay curve (time is not very small, not large) we assume 

F(T) = exp(- 7 T) = 1 - 7 T, 7 = const. (8) 
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Then the quantity T = 2Re7 is the probability of decay per unit of time (decay constant). 
Using Eq. (^) and Eq. (§), we obtain 

roc 

T = 2Re / (^ \V'e- i ^ +w ^V'\^ )e i£ o t dt. (9) 
Jo 

By v(r), Ey) denote the matrix elements of V which cause the decay of state \X e ) and emitting 
of particle Y: 

v(rj,E Y ) = {X 9 Y(r 1 ,Ey)\V X y\X e Y ). (10) 
All other matrix elements don't effect on the decay constant. Let us introduce the vector 

\Y) = J d V dE Y \Y(ri, E y ))v(t), E y ). (11) 

Then, after simple algebraic transformations, Eq. @ may be rewritten as 

r = 2tt J™ M(E Y )A (Ey - E Y in ) dE Y , (12) 

where 

E Y m = Ey+u + 5V , u Q = E e x -E x , 



M(E Y ) = Jdr)\v(r),Ey) 



|2 



A(E) = -Re / " D{t)e~ iEt dt, (13) 



7T 



(YZ 




{H°, z +W YZ )t 


YZ ) 


(YZ 




YZ ) 



Hy Z = Hy®I Z + Hl®Iy. 

It is easily shown that / A(E)dE = 1. We can not obtain an explicit analytical expression 
for r with respect to the matrix elements of interaction W due to nonperturbative character 
of this interaction. However, it is possible to derive an interesting qualitative conclusion on 
the shape of function A(E) (which is essentially used to calculate T) without calculations. 
For simplicity we suppose M(Ey) to be a step-like function with the jump at energy E Y r . 
Suppose the energy-spreading function A(E) be a bell-like and nearly symmetric with the 
maximum at E = 0. Consider the case E Y in < Ey 11 * (Fig.D. Then the transition \X e ) -> \X a ) 
with emitting of particle Y is strictly forbidden by the Energy Conservation Law. For 
example, the binding energy of atomic electron is greater than the transition energy uq and 
the electron can not be ionized during this transition. But the functions A(Ey — E Y m ) and 
M(Ey) may have no- zero overlap integral Eq. (|1^) as is shown on Fig. 0. Hence T > and the 
transition \X e ) — > \X g ) is possible. Thus, we come to a contradiction. This contradiction 
means that the suggestion about shape of function A(E) was wrong. The contradiction 
could be eliminated if whole no-zero part of function A(E) locates at the left-hand side of 
point E = 0. Therefore, we conclude, that our formalisms predict this special shape for 
the function A(E). Another prediction is the following. Let E Y m — Ey r > but the value 
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Figure 2: No-zero probability of decay forbidden by the Energy Conservation Law. 



Ey — Ey r is of the order of function A(E) width or less. Then the transition \X e ) — *> \X g ) 
is permitted, but a considerable part of function A(Ey — E Y m ) is located at the left-hand 
side of point Ey r , so T < T . Here T is the decay constant not perturbed by the interaction 
of particle Y with the device Z. This is QZE. In the following sections we will verify both 
predictions by direct calculations with a simple numerical model. 



3 Numerical model 

We consider one-dimensional three-particle model (Fig.|3|) in this section and hereafter in 
present paper. The systems X, Y, Z are one- dimensional rectangular potential wells. There 
is a single particle in each well in the initial state of system X <g> Y <S> Z. The masses of 
particles and the geometry of potential wells are clear from Fig.|3|. We use the units such 
that my = 1, cly = 1, fr = 1. The coordinates of particles X, Y, Z are denoted by 
x, y, z, respectively. There is infinitely high potential wall for all particles at the point 
x = y = z = 0, consequently all particle eigenstates are no-degenerated. We consider that 
each particle X, Y, Z governs only by its own potential well Ux(x), Uy(y), Uz(z) respectively 
and by interparticle interactions. 

The potential well Ux{x) is a potential box with solid walls. The potential wells Uy{y) and 
Uz{z) are such that they contain only one bounded state for particles Y and Z, respectively. 
The particles X and Y interact by repulsive £-like potential 

Vxv(y - x) = v 8(y - x), v >0. (15) 

This interaction causes transition of particle X from the initial state \X e ) to the ground 
state \X g ) and simultaneously excitation of particle Y from the bounded state \Yq) to the 
continuum \Y b (Ey))- Since all states are no-degenerated, the degeneration index rj may be 
omitted. The threshold energy for particle Y to be ionized is Ey r = 0. The particles Y and 
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Figure 3: The one-dimensional three-particle model of spontaneous decay with continuous 
observation of decay particle by distant detector. 



Z interacts by Gaussian repulsive potential 

W Y z(z -y) = Pq w exp 



1< 



Pq, w > 0; 



(16) 



P^ = (I YZ -\Y Z )(Y Z \). 



The potential Wyz always fulfills the condition Eq. (fj) due to the artificial factors Pq. We 
discuss these factors in the last section of paper. The Hamiltonian of joint system X <g> Y ® Z 
is 



H 



1 d 2 



2m T dx 2 



+ U x (x 



Iyz + 



1 d 2 



2m y dy 2 



+ U Y {y) 



L XZ 



+ 



i o 2 



+ U z {z) 



® I X Y + V X y{y -x)®I z + W YZ {z -y)® Ix- 



2m z dz 2 

General expressions for v(t],Ey), \Y), and T (Eqs. (|10| ,pLT],pL2D respectively) now become 

v(E Y ) = (X g Y(E Y )\V XY \X e Y ), (17) 
\Y) = f \Y(E Y ))v(E Y )dE Y , (18) 



r = 2tt f \v{E Y )\ 2 A(E Y - Ep n )dE Y . 



(19) 



The expressions for A(E) and D(t) (Eqs. (|13],[l4|) respectively) are remain unchanged. 

To calculate T we should calculate D(t). To calculate D{t) we should calculate two 
functions 



(YZ \e-< H rz +WYZ »\YZ ), 



q (t) = (YZole-^lYZo) 



(20) 
(21) 



and then find D(t) = q{t)/qo(t). In this paper we calculate q(t) numerically. 
To calculate q(t) the Shrodinger equation may be solved: 



M(y,z,t) 
*(y,z,0) 

and then the inner product q(t) = 



(H YZ + W YZ )^(y,z,t) 
Y(y)Z (z) 



(22) 
(23) 



Eqs. ([ID 

Y as 



(YZo\ty(y, z, t)) may be obtained. It follows from 
(|l7l) , and (|j|) that Y(y) may be represented through functions X e , X g , and 



Y(y) = NY (y) 



X*(y)X e (y)- / \Y {y')\ 2 X:{y')X e {y')dy' 



(24) 



where N is a normalization factor. Since the functions X e (x), 



X g (x) 



Y (y), and Z (z) are 



well known eigenf unctions of one-dimensional rectangular well, it is easy to calculate the 
initial state Eq. fl23| ) analytically. Note that it follows from Eq. fl24"|) that Y(y) is a compact 
wave packet near the origin of axis y. The physical meaning of this wave packet is that it is 
the particle Y state that arises virtually just after the particle excitation |T2| . 

Eq. (|22|) was solved numerically. The state of the system Y <g> Z was represented by 
a grid wave function with zero margin conditions defined on two-dimensional equidistant 
rectangular grid with the same steps along y- and z-axis. Both dimensions Ly and Lz of 
calculation area were much greater than distance zq from the center of device Z to the origin 
of coordinate system. The scheme of calculation was as follows. Let the grid wave function 
at the time t be {^^(t)} where k = 0, . . . , Ny~, 1 = 0,..., Nz- Then the wave function at 
the time t + At is calculated through successive four steps (a), (6), (c), (d): 
(a) Calculation of sin-Fourier transform of the grid function {^^(t)}: 



N Y -1 N z -1 

E E 

k=l 1=1 



~ , . . (mix , 
sin ^— k 



N Y N Z 

(b) Calculation of free evolution of Fourier coefficients: 



sin 



F mn (t + At) = F mn (t) exp 



1 



/mil 



2my \ L 



Y 



2m z 



nit 



nix 



At 



(c) Calculation of back sin-Fourier transform that produces the free evolution of system 
Y ® Z without potentials Uy, Uz, and Wyz during time interval At: 



*«(* + Af) 



N Y -1 N z -1 

E E 

k=l 1=1 



Finn [t 



. . . frrm , 
Ai) m (_fc 



sm 



8 



Table 1: The parameters of problem for the models "Wide W" and "Narrow W" . Here U Y 
and U z are the depths of wells Uy and U z . 



Parameter 


Wide W 


Narrow W 


a x 


0.6 


0.6 


Uly 


1.0 


1.0 


ay 


1.0 


1.0 


Uy 


-5.552 


-5.552 


H,y 


-2.776 


-2.776 


zo 


4.0 


4.0 


2a z 


1.0 


1.0 


m z 


0.9 


0.9 


u° z 


-2.210 


-2.210 




-1.0 


-1.0 


a w 


2.548 


0.2 


Wq 


789.2 


20000 



(d) Calculation of contribution of all interactions to the evolution during time interval At: 

V kl (t + At) = *' w exp{-i[£7y(2/ fe ) + U z (zi) + W{z { - y k )]At}. 

The zero margin conditions is fulfilled because of representation of {^m} by sin-Fourier 
series. 

The calculation of function qo(t) Eq. ([?!]) is not difficult. This calculation may be carried 
out analytically or numerically by the same way as the calculation of function q(t) but for 
Wy Z = 0. To verify our calculation schemes both ways was tested (the results was identical). 

To calculate the function A(E) through D{t) one should calculate the Fourier transform 
Eq. flTSp. To do this we used the cubic spline approximation of the numerical function D{t). 

4 Results of calculations and discussion 

We present the results of calculations for two sets of parameters of problem (Table[IJ). All 
parameters were the same for both calculations except the parameters of interaction W. The 
first variant was the "Wide W" . For this variant W was wide enough for particle Z in its 
ground state \Z ) to feel the appearance of the particle Y in the continuum spectrum near 
the origin of coordinate system. Also, W was strong enough to ionize the particle Z from 
its ground state. The second variant was the "Narrow W" . For this case W was narrow 
enough that the particle Z does not feel the particle Y near the origin of coordinate system. 
Also, W was strong enough for particle Y could not be tunnelled through particle Z and the 
energy of transition was high enough to ionize Z. 

The results of "Wide W" calculation is presented on Fig. [|. One can see from Fig. 
that the function q(t) drops down faster then the function qo{t). This is a result of detector 
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Figure 4: The calculations for "Wide W" . (a): solid line — \q(t)\, dashed line — |<7o(£)|; (b): 
solid line — dashed line — ReD(t), dotted line — Im D(t); (c): A(E); (d): solid line- 
dependency of unperturbed value of decay constant T on final energy Ey ', dashed line- 
dependency of perturbed value T on Ey. T(E Y m ) ~ for < Ey < 100 (Zeno effect). 



Z excitation due to the interaction Wyz between the particle Y and the detector Z (see 
Eqs. (|20"1J2ID ). As a result, the absolute value of function D{t) = qo(t)/q(t) drops down from 
the initial value 1.0 to zero (Fig. At the same time the imaginary and real parts of D(t) 
oscillate. The function A(E) (Fig.|]c) is a Fourier transform of D(t) (Eq. (|I3|) ), therefore this 
function is bell-like due to dropping of function D(t) and has left-side shift due to oscillations 
of real and imaginary part of D(t). Moreover, it is seen from Fig. that A(E) « for E > 0, 
as it was predicted on the base of Energy Conservation Law in Section 0. 

One can change the transition energy ujq of system X by altering the particle mass nix- 
Then the final energy Ey w uq — \E Y \ of particle Y is changed simultaneously. Therefore, 
it is possible to consider the dependency of decay constants on Ey . The dependency of 
unperturbed decay constant r (i. e. for Wyz = 0) on E Y m is shown on Fig.^rf by solid 
line. The complicated shape of this function is a consequence of particle Y reflection from 
the sharp margins of Uy potential well. The dependency of decay constant F perturbed 
by interaction Wyz on Ey is shown on Fig.£|(i by dashed line. It is seen that T{E^ in ) is 
strongly inhibited in comparison with r for values of E Y m which are less then approximately 
100. This is the predicted in Section |2| Zeno effect. Zeno effect in spontaneous decay takes 
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Figure 5: The calculations for "Narrow W" . Solid line — D(t); dashed line — P sur (t). 

place for low energies of decay particles, near the threshold of decay, if this effect is presented 
at all. 

The model with "Wide W" interaction does not contradict any fundamental principles 
of quantum theory, but this model is quite unrealistic practically. The long distance inter- 
action between Y and Z must couple the ground states \Y ) and \Z ) of systems Y and Z 
inevitably. To obtain Wyz\YoZ ) = we inserted the artificial factors [lyz — \YoZo)(Y Z \) 
into the interaction Wyz i n Eq. (|16D. It would be more realistic to consider sufficiently 
narrow interaction Wyz to obtain by natural way 



w exp 



yf 



Y (y)Z (z) « 0. 



Then the factors (lyz — \YoZo) (Y Zo\) may be omitted, they do not play a role any more. 
To consider this realistic situation we studied the model of "Narrow W" with aw = 0.2. It 
is seen that aw *C Zq. The results of calculation with "Narrow W" was quite different from 
"Wide W" ones. The function q(t) occurred to be almost the same as function qo(t). The 
resulting function D(t) is shown on Fig. [5] by solid line. It is seen that D(t) does not show 
a drop-down behavior, but rather shows some oscillations at long times. It is impossible 
to calculate the function A(E) numerically in this situation, because the integral Eq. (|l"3|) 
diverges, but it is clear that A(E) will be 5-like, not spreaded bell-like function. Thus, T 
and To are almost equal each other and Zeno effect is absent in "Narrow W" model. 

We mentioned that it would be reasonably to consider the function A(Ey — E Y m ) in 
Eq. ([12]) as an energy spreading of final state of decay due to a finite time life of particle 
Y until inelastic scattering on detector Z. Then the sense of function D(t) is the effective 
"decay curve" of the final state of decay in the analogy with the decay onto an unstable level 
|T3| . But it is clearly seen that it is not the case for the model of "Narrow W" . Survival 



probability of the state \Z ) after decay of the system X was occur, may be written as 



P sur (t) = Tr [\Z )(Z \p z (t)] = J dy 



dz^(y,z,t)Z*(z) 
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where ^{y,z,t) is the solution of Eq. (p0|) and pz(t) is the reduced density matrix of the 
system Z. The curve -P sur (£) is shown on Fig. [5] by dashed line. It is seen that the survival 
probability decreases with time (as could be expected) and that P SU r{t) is quite different 
from the function D{t). 

A general cause of the found D(t) behavior is the following. Note that for the model of 
"Wide W" the cause of D(t) dropping down is excitation of system Z. We would expect 
that the same reason may lead to dropping of D(t) in the model of "Narrow W" as well. 
Now consider the right hand side (RHS) term in Eq. (|20|). The function Y(y) is a compact 
wave packet near the origin of coordinate system. This packet contains both low-energy and 
high-energy components. Hence, the wave packet Y(y) does not drive with time from the 
left to the right along the y-axis, but spreaded out by the manner that the lowest energy 
part retains near the origin of coordinate system forever. This is a consequence of my > 0. 
Namely this lowest energy part of wave packet determines the value of inner product in 
the RHS of Eq. (|20| ) and, consequently, behavior of D(t). But this lowest energy part of 
wave packet can not influence on the state of system Z due to a short range of interaction 
Wyz(z — y) and large length of distance zq (Fig. ^). As a result there is no influence of system 
Z excitation to the behavior of function D{t). Oscillations of D{t) (Fig. [5]) is appeared to be 
a result of elastic reflection of particle Y from the particle Z in its ground state \Zq). 

Thus, our conclusions are as follows. Firstly, the functions A(E) and D(t) in Eqs. ( |12| . |l"3| ) 
have no any simple physical sense in the context of problem of continuous observation of decay 
by distant detector. Generally, the function D(t) does not mean the survival of the final state 
of decay with time generally, and the function A(E) does not mean the energy spreading 
of decay final states. Secondly, there is no quantum Zeno effect during the continuous 
observation of spontaneous decay by distant detector if interaction between emitted particle 
and detector is short-range and the emitted particle has no-zero mass. Thirdly, Zeno effect 
in the context of the same problem takes place if the interaction between emitted particle 
and detector is long-range, but this situation is considered as unrealistic. Finally, we did not 
consider the case of "spreaded" detector, when the detector is represented by some medium 
which contains a decay system and we did not consider the case of massless emitted particles. 
The existence of Zeno effect in these situations is meanwhile an open question. 
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